Complete population transfer in a three-state quantum system by a train of pairs of 

coincident pulses 
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A technique for complete population transfer between the two end states |1) and |3) of a three- 
state quantum system with a train of N pairs of resonant and coincident pump and Stokes pulses is 
introduced. A simple analytic formula is derived for the ratios of the pulse amplitudes in each pair 
for which the maximum transient population Pzit) of the middle state \2) is minimized, P™ ax = 
sin 2 (7r/4iV). It is remarkable that, even though the pulses are on exact resonance, Piit) is damped 
to negligibly small values even for a small number of pulse pairs. The population dynamics resembles 
generalized 7r-pulses for small N and stimulated Raman adiabatic passage for large N and therefore 
this technique can be viewed as a bridge between these well-known techniques. 
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Introduction. Coherent excitation of a discrete quan- 
tum system by an external resonant field represents 
an important notion in quantum mechanics. Resonant 
pulses of specific pulse areas are widely used in a vari- 
ety of fields, including coherent atomic excitation [H, 2|, 
nuclear magnetic resonance Q , quantum information 4 j , 
and others. Resonant excitation allows one to establish 
a complete control over the quantum system, particu- 
larly in two- and three-state systems. Important exam- 
ples include complete population transfer between the 
two states in a two-state system and between the two 
end states in a three-state A-system @, Q . Crucial con- 
ditions for resonant excitation are exact pulse areas and 
exact resonances between the frequencies of the external 
fields and the Bohr transition frequencies. Deviations 
from exact resonances or exact pulse areas lead to devia- 
tions of the transition probability from the desired value. 

To this end, an alternative to resonant excitation is 
provided by adiabatic passage techniques, which are ro- 
bust against such deviations. In three-state A-systems, 
the famous technique of stimulated Raman adiabatic pas- 
sage (STIRAP) allows complete population trans- 
fer between the two end states |1) and |3) in the adia- 
batic limit without placing any transient population in 
the middle state |2), even though the two driving fields 
- pump and Stokes — can be on exact resonance with 
their respective transitions, |1) -0- |2) and |2) O |3). The 
conditions for STIRAP are the two-photon resonance be- 
tween states |1) and |3), the counterintuitive pulse order 
(Stokes before pump), and adiabatic evolution. However, 
adiabatic evolution requires large pulse areas, typically 
over 107T for smooth pulse shapes, which may be hard to 
reach experimentally. Strategies for optimization of STI- 
RAP with minimal pulse areas have been developed 
but these come on the expense of strict relations on the 
pulse shapes [8|, [9( and require specific time-dependent 
nonzero detunings Q. 

In the present paper we introduce a novel technique 
that is both an alternative to the above techniques and 
reduces to them in two opposite limits. The technique en- 
ables complete population transfer between the two end 



states |1) and |3) of a three-state A-system with a train 
of N pairs of resonant and coincident pump and Stokes 
pulses with negligibly small transient population in the 
middle state |2), which vanishes as ./V -2 with the number 
of pulse pairs N. We note that unlike the adiabatic solu- 
tions, which are approximate, our technique is described 
by an exact analytic solution. This technique formally 
resembles the techniques of piecewise adiabatic passage 
[ill [l2| and composite pulses [H, [H| ; the differences will 
be discussed toward the end. 

We shall first present the exact analytic solution to the 
three-state dynamics for a single pulse pair and then the 
exact solution for a train of N pairs of pulses, in which we 
shall demonstrate explicitly the dynamical suppression of 
the middle-state population. 

Single pulse pair. The probability amplitudes Ck(t) 
(k = 1,2,3) of the three states that form the A-system 
obey the Schrodinger equation, 



ihd t c(t) = H(t)c(t), 



(1) 



where c(t) = [ci(t), C2(t), c^{t)] T . The Hamiltonian in 
the rotating- wave approximation [l], Q reads 



H(t) 
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where fl p (t) and £l s (t) are the Rabi frequencies of the 
pump and Stokes pulses, respectively; each of them is 
proportional to the electric-field amplitude of the respec- 
tive laser field and the corresponding transition dipole 
moment, fl p (t) — — di 2 • E p (t) and Sl s (t) = — d 32 • E s (t). 
r is the rate of irreversible loss from middle state |2). For 
simplicity both Q p (t) and Q. s (t) will be assumed real and 
positive because their phases can be eliminated by redef- 
inition of the probability amplitudes. More importantly, 
we assume that the Rabi frequencies are pulse-shaped 
functions with the same time dependence f(t), but pos- 
sibly with different magnitudes, 
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In this case — single- and two-photon resonances and the 
same time dependence of the pump and Stokes fields — 
the Schrodinger equation ([T]) is solved exactly by making 



a transformation to the so-called bright-dark basis [1 
The exact propagator reads [l(| 
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where the root-mean-square (rms) pulse area A and the 
mixing angle 9 are defined as 



A = 



tan 9 = 



n,(t) ~ ^ 



(5) 
(6) 



with A p s = J^ f Q P:S (t)dt being the pump and Stokes 
pulse areas. (For reasons that will become clear below 
we have omitted the area A from the arguments of U.) 
Due to condition ([3]) the angle 9 is constant. 

The propagator (0| allows us to find the exact ana- 
lytic solution for any initial condition; however, we re- 
strict our attention here to a system initially in state |1): 
c(ti) — (1, 0, 0) T . Then the populations at the end of the 
interaction are 



Pl(tf) 
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P 2 (tf) = sin 2 9 sin 2 
P 3 (t f ) = sin 2 29 sin 
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Obviously, complete population transfer to state |3) is 
achieved for 9 = n/4, which corresponds to A p = A s , and 
rms pulse area A — 2-k. Then, however, the intermediate 
state |2) acquires a transient population which reaches a 
maximum value of p™ ax = I at the intermediate time 
when the accumulated rms pulse area A(t) reaches the 
mid-point value it. This scenario is illustrated in Fig. [T] 
(left frames). 

We shall show below that the application of a train 
of resonant pulse pairs allows one to transfer the popula- 
tion from state 1 1} to state |3) completely, while reducing 
the transient population in state |2) to arbitrarily small 
value. 

Train of pulse pairs. A sequence of N pairs of pulses, 
each with rms pulse area A and mixing angles 9k, pro- 
duces the following total evolution matrix 



UW = U (9 N ) U (0 N -i) • • • U (9 k ) ■ • • U (0 X ) 



(8) 



When the system is initially in state |1), the final popu- 
lations are given by P n — \U^\ 2 , with n = 1,2,3. Our 



objective is to have P\ = P^ = and P3 = 1 at the end 
of the pulse train, with as little transient population in 
the middle state |2) as possible. 

It is convenient to have an "anagram" pulse train that 
is symmetric with respect to time reversal, i.e., with mix- 
ing angles 9 N+1 ^ k = tt/2 - 9 k , (k = 1, 2, 3... [N/2\ ). In 
order to determine the values of 9 k and A, we use Eqs. ^ 
and ((5J and we first demand that after each pulse pair 
the population of state |2) vanishes; this gives immedi- 
ately rms pulse area A = 2tt for each pulse pair, hence 
the omission of A from the arguments of U in Eq. (0| . 

Next we require that the population is transferred to 
state |3) in the end, Pi = P2 = and P3 = 1. We further 
require that the maximum of the transient population 
P2(t) excited by each pulse pair is the same. Among 
the many solutions we pick the one that minimizes the 
maxima of P2 (t) . A simple algebra gives the angles 



(2k- 1)tt 
AN 



(k = 1,2, 3. ..AO- 



(9) 



The maximum population of state |2), which, as for 
a single pulse pair above, occurs in the middle of each 
pulse pair (at rms area 7r), is readily obtained, 



pmax 
r 2 
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(10) 



From here we conclude immediately that for large N 
the maximum population of the middle state vanishes 
as 1/N 2 . Obviously, for N ^ 8 pairs, the transient pop- 
ulation in state |2) does not exceed 1%. It is particularly 
significant that this suppression occurs on resonance and 
it results from the destructive interference of the succes- 
sive interaction steps, rather than from a large detuning. 
We note that for N ^> 1 the total pulse area is very 
large, which is the condition for adiabatic evolution on 
resonance in STIRAP. 

We show in Fig. [T] the population evolution for sev- 
eral pulse trains of different number of pulse pairs N. In 
all cases the population is transferred from state |1) to 
state 1 3) in the end in a stepwise manner. The transient 
population of the intermediate state |2) is damped as N 
increases: from 0.5 for a single pair of pulses and 0.15 for 
two pairs to below 1% for 8 pulse pairs. For small N the 
transition picture resembles (fractional) resonant excita- 
tion, whereas for large N the transition picture resem- 
bles adiabatic passage. The lower frames demonstrate 
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FIG. 1: Rabi frequencies (top frames) and populations (bottom frames) vs time for (from left to right) 1, 2, 3, 5 and 8 pairs 
of pulses. The pulse shapes are Gaussian, 0, p (t) — Q,s'm8ke~ < • t ~ Tk,| / T and fi s (t) = f2 cos # fe e~( t_Tfc ) / T , with O, = l^pK )T 
(corresponding to rms pulse area A = 27r) and the mixing angles 9k (k = 1, 2, . . . , N) are given by Eq. (|9|. Middle frames: no 
decay (T = 0). Complete (stepwise) population transfer 1 — > 3 is achieved in all cases; however, the population P2(t) of the 
intermediate state is different. The maximum of Pzit) is given by Eq. (I10[l ; it vanishes as 1/N 2 . Bottom frames: irreversible 
loss from state \2) with a rate F = 1/T. 



the effect of irreversible population loss from state \2) for 
a loss rate T = 1/T. As it can be expected from the 
middle frames of lossless interaction, the damping of the 
intermediate-state population by longer pulse trains re- 
duces the effect of population loss: the target-state popu- 
lation P3 decreases to 0.63 for a single pulse pair but just 
to 0.93 for a train of N = 8 pulse pairs. This population 
loss can be decreased even further by longer trains. 

We note that because the technique uses fields on exact 
resonance the pulse shapes are unimportant. Although 
the example in Fig. [T] uses Gaussian shapes, pulses of 
rectangular or any other shape are equally suitable. 

Discussion. We have demonstrated that complete 
population transfer in a three-state system driven by a 
pair of external pulse-shaped laser fields that are on reso- 
nance and have the same time dependence, can be accom- 
plished in such a way that all population is transferred 
from the initial state |1) to the target state |3) with min- 
imal transient population in the middle state |2). This is 
achieved with a train of TV pairs of coincident pulses with 



appropriately chosen amplitudes. In the limit N 3> 1, the 
present technique resembles STIRAP, because then the 
successive increments of the mixing angle 6k become very 
small (nearly continuous), with the Stokes field dominat- 
ing over the pump field in the beginning and then the 
pump field dominating in the end, in exact analogy to 
the counterintuitive sequence Stokes-pump in STIRAP 
0-0] • However, the present technique achieves complete 
population transfer 1 — > 3 also for small N, in a man- 
ner reminiscent of generalized 7r-pulses [2| , while keeping 
the middle-state population at very low values. In this 
manner, the present technique can be viewed as a bridge 
between generalized 7r-pulses and STIRAP. 

The present technique can be viewed also as an alterna- 
tive of two other techniques that use pulse trains for com- 
plete population transfer. Piccewise adiabatic passage 
(PAP) [ll], OjJ uses a train of a large number of pulses, 
each of which produces a perturbatively small change in 
the populations, while the present technique works for an 
arbitrary number of pulse pairs and each pair may pro- 
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duce a large population change (for small N). In one of 
the implementations involving only two states [l2[ , PAP 
demands phases that change quadratically from pulse to 
pulse, which translate into a linear chirp for a large num- 
ber of pulses; the present technique does not need such 
quadratic phases but appropriate amplitude ratios. In 
the PAP implementation with three states [ll| , to which 
the present technique is more closely related, the pump 
and Stokes fields in STIRAP are turned abruptly on and 
off repeatedly; the amplitudes of the individual pulses 
are determined such that they match the envelopes of 
the pump and Stokes pulses in STIRAP. In the present 
technique the pulse amplitudes are determined from the 
conditions to achieve complete population transfer to the 
target state |3) and to minimize the population of the 
middle state |2). Indeed, the systematic suppression of 
the transient population of state |2) with the increasing 
number of pulse pairs N seen in Fig.[T]is only observed in 
the present technique. We also point out that the solu- 
tion in the present paper is exact while PAP and STIRAP 
give approximate solutions in the adiabatic limit. 

The present technique is also reminiscent of the tech- 
nique of composite pulse sequences [l~3T ] . The latter uses 
sequences of pulses in two-state systems 1131. or sequences 
of pulse pairs in three-state systems [lj], with well- 
defined relative phases, which are determined from the 
condition to produce a desired excitation profile. There- 
fore the control parameters in the composite pulses are 



the relative phases, whereas in the present technique the 
control parameters are the amplitude ratios in each pulse 
pair. Moreover, the objective in the composite pulse 
technique is the shape of the excitation profile while in 
our technique the main objective, beside the complete 
population transfer 1 3, is the suppression of the 
intermediate-state population. 

Conclusion. The technique introduced in this pa- 
per allows complete population transfer between states 
|1) and |3) via an intermediate state |2) with a train of 
N pairs of coincident pump and Stokes pulses, by plac- 
ing only a negligible transient population in state |2), 
which decreases as 1/N 2 as the number of pulse pairs 
N increases. This technique resembles the technique of 
generalized 7r-pulses for small N and STIRAP for large 
N and therefore it can be viewed as a bridge between 
these two well-known techniques. It is remarkable that 
the middle-state population P2{t) is damped consider- 
ably even for a small number of pulse pairs despite the 
fact that the pump and Stokes fields are on exact res- 
onance with their transitions. All these features make 
this technique an interesting alternative of the existing 
techniques for coherent control of three-state quantum 
systems. 
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